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Abstract 
We obtain the reduced CMB data (l4, R, z.] from WMAP9, WMAP9+BKP, Planck+WP and 
Planck+WP+BKP for the ACDM and wCDM models with or without spatial curvature. We 
then use these reduced CMB data in combination with low-redshift observations to put constraints 
on cosmological parameters. We find that including BKP results in a higher value of the Hubble 
constant especially when the equation of state of dark energy and curvature are allowed to vary. For 
the ACDM model with curvature, the estimate of the Hubble constant with Planck+WP-+Lensing 


is inconsistent with the one derived from Planck+WP-+BKP at about 1.2 ø confidence level. 
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I. INTRODUCTION 


Since the discovery of the cosmic acceleration expansion of the universe based on the 
distance measurement of type Ia Supernovae (SNe Ia) Al d. its origin has become a hot 
topic in modern cosmology and theoretical physics. The cause to the observed cosmic 
acceleration is due to the so-called dark energy with negative pressure in general relativity 
framework, or the modification to general relativity at cosmic scales. 

To study properties of dark energy, one may combine some mature probes, such as SNe Ia, 
the observational Hubble parameter (HUB), baryon acoustic oscillation (BAO) and cosmic 
microwave background (CMB) anisotropy. The SNe Ia, HUB and BAO probe the expansion 
of the universe at low and intermediate redshifts, while the CMB measurements probe the 
distance at high-redshift (especially the distance to the surface of last-scattering). The CMB 
data provide the strongest constraints on cosmological parameters |3| and help break the 
degeneracies among the dark energy and other cosmological parameters. 

The reduced CMB data (l4, R, z.) provide an efficient summary of CMB information, 
of which /4 is the angular scale of the sound horizon at recombination and determines the 
acoustic peak structure of the CMB angular power spectra, R is the scaled distance to the 
recombination and determines the amplitude of the acoustic peak and z, is the redshift at 
the last scattering surface. Instead of the full CMB spectra, the reduced CMB data relating 
the distance to the last scattering surface provide a fast and self-consistent approach for 
combining the CMB information with complementary cosmological data to constrain late- 
time cosmological parameters. The reduced CMB data were firstly derived in |A| from the 
three-year WMAP data for the ACDM model (with and without spatial curvature) and 
they found that dark energy density is consistent with a constant in cosmic time and a flat 
universe is allowed by using the reduced CMB data together with the SNe Ia and BAO 
data. Recently Wang and Wang obtained the reduced CMB data from WMAP9 and 
Planck4-WP--Lensing for the Q,+ACDM model and found that the reduced CMB data 
derived from Planck+WP-+Lensing are much tighter than those from WMAP9, but when 
combined with other low-redshift observational data, the reduced CMB data do not improve 
the constraints on the dark energy too much compared to those from WMAP9. Besides, 
Shafer and Huterer ld derived the reduced CMB data from WMAP9 and Planck-- WP 
respectively for the flat wCDM model, when combining the reduced CMB data, BAO and 3 
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samples of SNe Ia data (Union2.1 "n SNLS3 [8| and PS1 "m respectively, they found that 
there is a preference for the equation of state of dark energy w « —1 for the constraints with 
Planck but not with WMAP9. More recently, the reduced CMB data have been obtained 
based on the flat (» CDM model by use of the Planck 2015 temperature and low-/ polarization 
data, showing that the reduced CMB data are consistent with those based on ACDM and 
CPL models " There are many works that use the reduced CMB data together with 
complementary cosmological data to constrain late-time cosmology 14. 

There are two advantages of the reduced CMB data. First, the CMB shift parameters 
[4 and R together with the decoupling redshift z,, extracted from the CMB angular power 
spectra, allow one to quickly evaluate the likelihood of various dark energy models, without 
the need to run a Markov Chain Monte Carlo exploration of the CMB likelihood which 
usually includes a number of astrophysical parameters to describe unresolved foreground 
components and other nuisance parameters. Second, it provides an efficient and appropriate 
summary of CMB data as far as dark energy constraints are concerned. Since l4 determines 
the acoustic structure in CMB angular power spectra while R determines the overall am- 
plitude of the acoustic peaks, they are nearly uncorrelated. Both R and /4 can be used to 
further compress CMB information and combined with other measurements in a friendly 
user manner to constrain dark energy models. 

On the other hand, the authors of analyzed the likelihood of the reduced CMB 
data with WMAP3 data for the base ACDM model involving extra parameters, such as 
tensor modes and a running spectral index. They found that adding curvature or slightly 
modifying the dark energy parameters does not significantly change the values of (l4, R}, 
which, however, change large when more parameters like tensor modes or running of the 
scalar spectral index are involved. 

The purposes of this work are to update constraints on the parameter combination 
UA, R, z4} using newly released CMB temperature and polarization data for several cosmo- 
logical models and to test their dependence on model assumptions. These updated reduced 
CMB data provide a simple and efficient method for combining in a friendly user manner 
the current CMB measurements with low-redshift data. We first obtain the reduced CMB 
data from WMAP9 data and Planck data together with WMAP polarization data 
(WP), based on the ACDM model and wCDM model with a flat or curved space curvature, 
respectively. We also make a joint analysis of the data from BICEP2/Keck Array |18] and 
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Planck (BKP) to derive the reduced CMB data. Our goal is to see the differences among 
the data used when combining the reduced CMB data with the low-redshift observational 
data to constrain different cosmological models. 

The paper is organized as follows. In section [I] we present the reduced CMB data ob- 
tained from WMAP9, Planck+WP, Planck4- WP--BKP and WMAP9-4-BKP based on differ- 
ent cosmological models, respectively. In section [III] we give the results of the combination 
of reduced CMB data and other data sets to constrain the different cosmological models. 


The conclusions are included in section 


II. REDUCED CMB DATA 


The distance measurement is one of the most powerful methods to study the evolution 
history of the universe. In a Friedmann-Robertson-Walker universe, the comoving distance 


from an observer to redshift z is given by 
r(z) = Hg* Q4 "2 sinn[|Q,| ^? T(2)], 


z dz! 
ro [zu EGGS. 


(1) 


where (2, = —k/H@ (k is the spatial curvature constant) and sinn(x) = sin(x), x, sinh(z) for 
Oy < 0, Qg = 0, and Qg > 0, respectively. The Hubble parameter is given by the Friedmann 


equation 


H?(z) = Hg [Q,o(1 + z)*+Qamo(1 + 2)? + woll + z^ 4 (2) 
O,(1 + z)? + (1 — Og — Qro — 04)], 


for the ACDM model, where the redshift z is defined by (1 + z) = 1/a, and 2,0, Qamo 
and Qo are the present values of the fraction energy density for radiation, dark matter and 
baryon matter, respectively. The latter two are often written as the total matter ld. 

16]: 


Omo = Qot Namo. The radiation density is the sum of photons and relativistic neutrinos 
Qro = QP (1 + 0.2271Nef f), (3) 


where Neff = 3.046 is the effective number of neutrino species in the Standard Model of 


particle physics bd. and QP = 2.469 x 10-5h-? for Temp = 2.725K (h = Ho/100 km s7! 
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Mpc 1). For the wCDM model, the Hubble parameter is given by 


H?(z) = HB[Q,0(1 + z)*+Qamo(1 + z)? + Qw(1 + 2)? 
Qx(1 + z)? + (1 — Ano — Qro — OK) F (z)], 


(4) 


where the evolving function F(z), depending on the equation of state of dark energy, is given 
by 

F(z) = (1+ 2*9. (5) 

It is noticed that here we have not assumed a flat universe model. As for a flat universe, 

the curvature terms disappear in equation (2) and equation (4) since Q; = 0. There are 

four cases considered here, which are the ACDM model and the wCDM model with a flat 


and curved space curvature, respectively. 


Data flat ACDM]|Q;,+ACDM flat wCDM |Q;--wCDM 
PLANCK+WP this work |Ref, [5] (--Lensing) Ref. [6] this work 


Table 1: References for the reduced CMB data derived from different CMB data in different 


cosmological models. 


In the CMB measurement, the distance to the last scattering surface can be accurately 
determined from the locations of peaks and troughs of acoustic oscillations. There are two 


quantities: one is the “acoustic scale" 


TD A(2,) 


la = (1+ 2,) AE 


; (6) 


and the other is the "shift parameter" 


R= 4/Qmo HE (1 + 2,)DA(z,). (7) 


Here D4(z) = r(z) bl +z) is the angular diameter distance and z, is the redshift at the last 
1 


scattering surface |2 
Z = 1048[1 + 0.00124(Qro A N OPIL + gi (Qmoh?)”I, (8) 
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Data + J J Correlation Matrix 


1.0000 0.3859 


WMAP9 - 0.66|1.7257 + 0.0165 |1088.96 = 0.3859 1.0000 


0.4998 0.8432 


1.0000 0.5126 
PLANCK-4-WP|301.66 + 0.18|1.7500 + 0.0089 |1090.33 4 0.5126 1.0000 


0.4552 0.8699 


1.0000 0.5526 


PLANCK+WP 
- 0.18|1.74974 + 0.0087 | 1090.04 + 0. 0.5523 1.0000 
+BKP 
0.4851 0.8725 
1.0000 0.4112 
WMAP9+BKP]301.83 + 0.66| 1.7210 + 0.0165 |1088.61 + 0. 0.4112 1.0000 


0.5217 0.8550 


Table 2: The mean values, standard deviations of {l4, R, z.) and the correlation matrix for the 


flat ACDM model. 


where 
7 0.0783 (Qroh?) 70238 


N= 14 39.5(Qoh2)9765" (9 
o 0.560 
927 TF LOk AT 


These quantities can be used to constrain some cosmological parameters without need to 
use the full likelihoods of WMAP9 or Planck data K 

Based on the original idea proposed in |4], Hinshaw et al. obtained constraints on the 
parameter combination (l4, R, z,) from WMAP9 data based on the wCDM model without 
assuming a flat universe. Wang and Wang obtained the mean values and normalized 
covariance matrix of (14, R, Qygh?, ns} from WMAP9 and Planck+WP-+Lensing data, re- 
spectively, based on the Q,+ACDM model. Recently, Shafer and Huterer ld derived the 
related results about (l4, R, z.) from WMAP9 and Planck+WP data, respectively, based on 
the flat wCDM model. In this section, following Wang and Wang [5| we obtain the Markov 
chains i" the Markov Chain Monte Carlo sampler as implemented in the CosmoMC 

2 


package [22] and then derive constraints on the parameter combination {l,4, R, z4}. In our 
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Data Correlation Matrix 
1.0000 0.3883  —0.6089 
WMAP9 - 0.66|1.732'7 + 0.016410.02260 + 0.00053] | 0.3883 1.0000 —0.5239 
—0.6089 —0.5239 1.0000 
1.0000 0.5250 —0.4475 
PLANCK+WP 
t 0.18|1.740'7 + 0.0094:10.02228 + 0.00030] | 0.5250 1.0000 —0.6925 
+Lensing 
—0.4475 —0.6925 1.0000 
1.0000 0.5279 
PLANCK+WP 
t 0.19|1.7416 + 0.0097 1089.74 4 0.5681 0.8946 
+BKP 
0.5279 1.0000 
1.0000 0.5217 
WMAPO9-4-BKP 301.94 + 0.66] 1.7251 + 0.0167]1088.88 4 0.4112 0.8550 


0.5217 1.0000 


Table 3: The mean values, standard deviations of (/4, R, zx} and the correlation matrix for the 


Q,+ACDM model. 


analysis, we focus on four cosmological models listed in Table [| based on the six-parameter 


model in the case of the flat ACDM model, described by 
{Qyoh?, up, 9, T, A,, Nns}, 


where O, is the ratio of the sound horizon to the angular diameter distance at the pho- 
ton decoupling, 7 is the Thomson scattering optical depth due to reionization, A, is the 
amplitude of primordial curvature perturbations and m, is the scalar spectral index. The 
CMB data sets used in our analysis are listed in Table [| Here we emphasize that both the 
tensor-to-scalar ratio r and running of the scalar spectral index a, are allowed to vary if 
the BICEP2 B-mode polarization data are included, because the B-mode power spectrum 
from the BICEP2 experiment implies the detect of primordial gravitational wave at 7.0 c 
ignoring foreground dust [19 and allowing the running of the scalar spectral index recon- 
ciles the tension with the Planck constraints on r |23]. These results have been confirmed 
by other data on the same field from the successor experiment Keck Array [18]. However, 


it is argued in lod that given the uncertainties of the amplitude of the dust polarization 


T 
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Data + J J Correlation Matrix 


1.0000 0.4077 0.5132 
WMAP9 - 0.66|1.7302 + 0.0169 1089.09 4 0.4077 1.0000 0.8580 


0.5132 0.8580 1.0000 


1.0000 


PLANCK4-WP |301.65 + 0.18|1.7499 + 0.0088 |1090.41 4 0.5262 
0.4708 
1.0000 
PLANCK+WP 
- 0.18|1.7495 + 0.0087 | 1090.31 + 0. 0.5379 
+BKP 
0.4782 
1.0000 0.4293 
WMAP9+BKP 301.88 + 0.66|1.7227 + 0.0174|1088.69 + 0. 0.4293 1.0000 


0.5267 0.8722 


Table 4: The mean values, standard deviations of {l4, R, zx} and the correlation matrix for the 


flat wCDM model. 


at the BICEP2 frequency of 150 GHz one cannot say conclusively at present whether the 
B-modes detected by BICEP2 are due to gravitational waves or just polarized dust. By 
using genus statistics, the authors of Ref " claim to find the evidence for the primodal 
gravitational wave signal with r = 0.11 + 0.04. Planck team |27] released the polarization 
data from 100 to 353 GHz, extrapolation of the Planck 353 GHz data to 150 GHz gives 
a dust power, which is the same magnitude as reported by BICEP2. Recently, Ref. 
performed a joint analysis of BICEP2/Keck and Planck data and obtained an upper limit 
r « 0.12 at 9596 confidence, showing little evidence of detection of primordial gravitational 
wave. Therefore the tensor-to-scalar ratio is allowed to vary if the BKP data are included 
in our analysis. 

The mean values, the standard deviations and their correlation matrix of (l4, R, z,] 
(or (14, R, Qwh?}) for four cosmological models by using different data are summarized in 
Table ]to Table B| respectively. It is noticed that Wang and Wang [5] used Q;9A? instead of 
z,, which gives identical constraints by replacing Qh? with z,. Different from other cases 


in the third row of Table [| they also used Planck data together with Planck lensing. 
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Data + J J Correlation Matrix 


1.0000 0.4262 0.5391 
WMAP9 - 0.67|1.7246 + 0.0183 | 1090.88 4 0.4262 1.0000 0.8643 


0.5391 0.8643 1.0000 


PLANCK+WP 301.60 + 0.18|1.7442 + 0.0093 |1089.86 4 


PLANCK+WP 


- 0.19|1.7407 4 


+BKP 


WMAP9+BKP|301.97 + 0.65|1.7251 + 0.0173 |1088.89 4 


Table 5: The mean values, standard deviations of {l4, R, zx} and the correlation matrix for the 


Q,+wCDM model. 


From TableB]to Table[5] we see that the Planck data give tighter constraints on (14, R, 24} 
than WMAP9 in the same cosmological model. Including the BKP data does not change 
the results significantly but the standard deviations seem to be a little larger, this is be- 
cause tensor perturbations have been considered when we use Planck4-WP--BKP to obtain 
constraints on (l4, R, z4}. We also notice that there is some tension between the WMAP9 
data and Planck data constraining on (la, R,z.}. For example, in Table 2] the constraints 
on R (z,) are inconsistent at about 1.5 ø (1.6 e) when using WMAP9 and Planck-- WP. 
The estimates of (I4, R, z4} are consistent with each other within 1 ø for the ACDM model 
when we use Planck4-WP and Planck+WP+BKP data. 

Moreover, the difference of values of the reduced CMB data derived from same data for 
different models is not significant. As stated in Ref. id. curvature or slightly modifying the 
dark energy parameters does not significantly change the values of (l4, R}. For example, 
from Table[2]and Table B] we find that the WMAPSO data give values of (14, R} = (301.95 + 
0.66, 1.7257 + 0.0165} ({302.02 + 0.66, 1.7327 + 0.01641) for the ACDM model without 


(with) spatial curvature which shows no significant difference. 
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III. COSMOLOGICAL PARAMETERS 


In the previous section we have derived the reduced CMB data from WMAPS, 
Planck-- WP, Planck+WP+BKP and WMAP9+BKP for the ACDM model and the CDM 
model with and without spatial curvature, respectively. In this section, we focus on con- 
straints on the cosmological parameters for the corresponding cosmological models from 
reduced CMB data in combination with the low-redshift observational data including the 
Union2.1 SNe Ia sample, Hubble parameter and BAO data, which are described in Ap- 
pendix. The best-fitted values of Qmo and A for the ACDM model, Qmo, h and w for the 
wCDM model and their 68% confidence level (CL) errors are given by using the Markov 
Chain Monte Carlo analysis in the multidimensional parameter space in a Bayesian frame- 
work. The results are summarized in Table [6] to Table [9] and their likelihoods are shown in 
Figure [I] to Figure [4] respectively. 

From Table [6] we see that in the context of the flat ACDM model, the combination of 
Planck data favors a relatively higher value of Qmo and a lower value of h compared to 
the combination of WMAP9 data. However, the reduced CMB data from Planck+WP do 
not lead to significantly improve the constraint on dark energy together with low-redshift 
observational data, compared to the reduced CMB data from WMAP9 even though the 
Planck measures all of the CMB distance parameters (l4, R, z,) more precisely, whose errors 
are 2 — 3 times smaller. This is because the Planck data appear to favor a higher value of 
Qmo and a lower value of Ho in the standard six-parameter ACDM model, which are in 
tension with the magnitude-redshift relation for SNe Ia and recent direct measurements of 
Ho . The constraints with BKP data suppress the value of Qmo and raise the value 
of h. 'The tendence seems to appear in all the cosmological models we are considering here. 
These estimates of Qmo and h are consistent with each other within 1 ø CL, but are in 
tension with the results derived by Planck i3. 

As we can see from Table [7] in the context of ACDM model with spatial curvature, the 
constraints with Planck+WP-+Lensing give h = 0.6880*00095, which is inconsistent with the 
value h = 0.6988*00090 derived from Planck4-WP--BKP at about 1.2 ¢ CL. The Planck team 
gives h = 0.6781 from TT+LowP+Lensing and h = 0.6731 from TT--LowP, and the value 
of h may be enhanced by Planck lensing data , so we can conclude that the constraints 


with Planck-- WP in the ACDM model with spatial curvature may give a lower value of h, 


10 
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which deviates more from h = 0.6988. Furthermore, the WMAP9 data favor a positive 0; 
(an open universe), but the Planck data give a negative €); (a closed universe). However, 
there is no evidence for any departure from a spatially flat geometry in these three cases. 

A cosmological constant has an equation of state (EOS) w — —1. If we release the EOS 
w of dark energy, the constraints with Planck+WP give w = —1.0507*00209, as shown in 
Table[8] which favors the phantom region at 1 c CL. The combination with BKP data gives 
a relatively higher value of w (w = —1.0489*0025:), while the constraints with WMAP9 
(WMAP9+BKP) give w = —1.0180*00225 (w = —1.0155*0053,.). All the three cases are 
consistent with the ACDM model. 

WMAPS, Planck+WP, Planck+WP+BKP and WMAP9+BKP all give negative €); (a 
closed universe) for the Qy+wCDM model, as shown in Table [9] but are consistent with a 
flat geometry within 1 ø CL. 

We also give the constraints for the Q} + (CDM model by using four sets of the reduced 
CMB data from Planck-- WP--BKP derived from flat ACDM, Q + ACDM, flat wCDM and 
Ox + wCDM, respectively, together with the low-redshift observational data. The purpose 
is to see whether the constraints on the cosmological parameters are sensitive to the choice 
of the reduced CMB data derived from different cosmological models. Their likelihoods 
are shown in Figure We can see that the likelihoods for the parameters h and Qmo are 
almost the same in the four cases. As for the parameters 2; and w, the likelihoods show 
no significant difference between the flat ACDM and flat wCDM cases, and between the 
Oy + ACDM and Qg + wCDM cases, respectively. However, it can be seen clearly that the 
reduced CMB data derived from the cosmological model with spacial curvature can give 
much better constraints on 9; than those derived from flat cosmological model. But the 


values of €), and w are still consistent with each other within 1 o CL in the four cases. 


Ll 
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Data Sno — Jp 

WMAP9 0.7022* 0095 
Planck+WP 0.697570 oon 
Planck4-WP--BKP |0.296770-0109 |9,6973 0.0079 
0.7039+9:0997 


Table 6: Constraints with 1 c errors on Qmo and A for the flat ACDM model from SNe Ia, HUB, 


WMAP9+BKP 


BAO and reduced CMB data. 


WMAP9 (994600115 | 0.6904+9:0112 | 0.0007+9:0051 
0.2986+9:012 |9.6880 00090 | —0.0006+9:9054 
Planckt+WP+BKP  |0.2963*5 015. [0.6988 00192 | —0.0001 0:009 


WMAP9+BKP 0.2904+9:0199 0.6000 2-015 | —0.0013+9:9066 


Table 7: Constraints with 1 ø errors on (3,459, h and Qg for the Q;,+ACDM model from SNe Ia , 


Planck+WP+Lensing 


HUB, BAO and reduced CMB data. 


WMAP9 0.287840 oia [0.7043 Oi | - 0180 0/06 
Planck+ WP [0.203671 0.704518: [1.050749 369 
Planck+WP-+BKP 0.29300 0106 0.7054 0011s | — 1.04897 0.0501 


WMAP9+BKP  .|0.2866*00115 |0.7049*20126 | 1.015510 028% 


Table 8: Constraints with 1 c errors on Qmo, h and w for the flat wCDM model from SNe Ia , 


HUB, BAO and reduced CMB data. 


IV. CONCLUSIONS 


The reduced CMB data provide an efficient summary of CMB information, and can be 
used to constrain cosmological parameters instead of the full CMB power spectra. We have 
obtained the reduced CMB data from WMAP9 data and Planck date based on the ACDM 


model and wCDM model with a flat or curved space curvature, respectively. We have 


12 
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WMAPS a [o voor | 1.052687 |-0.0001 48 018 


Planck+WP 0.2945* 00102 | 0.707820 Gr74 | 1.06147 0.0745 | — 0.00137 00090 
Planck-- WP--BKP [0.2940 00152 0.70957 0:07 | 1.0576--0:0636 | -0.0029* 00051 


WMAP9+BKP  .|0.2936*90171|0.7105* 0.0203 | 1.0437* 00951 |—0.00667 0.0100 


Table 9: Constraints with 1 c errors on Qmo, h, w and Qg for the Qy+wCDM model from SNe Ia, 
HUB, BAO and reduced CMB data. 


Fig. 1: Marginalized posterior distributions for h (right) and Qmo (left) of the flat ACDM model. 


also used the BKP data together with the WMAP9 and Planck data to derive the reduced 
CMB data. We have found that the Planck data give tighter constraints on (l4, R, z4} than 
WMAP9 in the same cosmological model. While including the BKP data, the standard 
deviations seem to be a little larger because of additional free parameters, the tensor-to- 
scalar ratio and running of the scalar spectral index. 

We have combined these reduced CMB data with low-redshift observational data to con- 
strain the cosmological parameters for the ACDM model and the «CDM model. The reduced 
CMB data from Planck+WP do not lead to significant improvement to the constraint on 
dark energy together with low-redshift observational data, compared to the reduced CMB 
data from WMAP9. Including BKP data results in a higher value of the Hubble constant 
especially when the equation of state of dark energy and curvature are allowed to vary. 
For the Q;+ACDM model, the constraint from Planck+WP-+Lensing in combination with 
low-redshift observations gives h = 0.6880*00099. which is inconsistent with the value of 


h = 0.6988* 0.0079 derived from Planck4-WP--BKP at about 1.2 0 CL. The constraint on w 


with Planck+WP gives w = —1.0507*00:02, favoring the phantom region at 1 ø CL, for a 
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0.0 
-0.000 -0.015 -0.010  -0.005 0.000 0.005 0.010 0.015 
ù 


Fig. 2: Marginalized posterior distributions for h (top left), Qmo (top right) and Qj, (bottom) of 
the Q,+ACDM model. 


flat wC DM model. The constraint on w with WMAPO9 or Planck+WP+BKP is consistent 
with w = —1. We have also found that the constraints on the cosmological parameter Qy 
and w are a little sensitive to the choice of the reduced CMB data derived from different 


cosmological models. 
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Fig. 3: Marginalized posterior distributions for h (top left), Qmo (top right) and w (bottom) of the 
flat wCDM model. 


Appendices 


A. "Type Ia Supernovae 
In this work, we take the Union2.1 compilation "n which contains 580 SNe Ia data over 
the redshift range 0.015 € z < 1.414. The chisquare is defined as 


i, y te) = inea (10) 


X = 
E i=l oS (%) i 


where u% (z) is the measured distance modulus from the data and u*”(z) is the theoretical 


distance modulus, defined as 
p^ (z) = Slogigdz + Ho, joo = 42.384 — 5log;g h. (11) 


The luminosity distance is 
dr(z) = (1+ z)r(z), (12) 


where r(z) is the comoving distance defined in d The nuisance parameter po can 
33| : 


be eliminated by expanding Y? with respect to po as 
xn = A + 2Buo + Cus, (13) 
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Fig. 4: Marginalized posterior distributions for h (top left), Qmo (top right), Qg (bottom left) and 


w (bottom right) of the Q,+wCDM model. 
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to the different models used to derive the reduced CMB data from Planck+WP+BKP. 
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where 


N obs 
A= y [u^ (zi; Ho -0)-u : (zl? 
i ENTE ? 
i=1 SNN 
N th obs 
u” (zi; Ho = 0) — p^" zi 
mes 25 E E : (14) 
i=1 SNA 
= 1 
g= 
2 agy (i) 
The x2, has a minimum as 
Yan = A- BC, (15) 


In this way the nuisance parameter uo is removed. This technique is equivalent to performing 
a uniform marginalization over Ho a) We will adopt X$y as the goodness of fitting instead 


of Xone 


B. Observational Hubble parameter (HUB) 


In this paper we use 19 observational Hubble data over the redshift range: 0.07 < z < 2.3, 
which contain 11 observational Hubble data obtained from the differential ages of passively 
evolving galaxies Lah and 8 H(z) data at eight different redshifts obtained from the 
differential spectroscropic evolution of early type galaxies as a function of redshift " The 


chisqure is defined as 
N 


T 2 Bald ee l (16) 


where H;,(z) and Hos(z) are the theoretical and observed values of Hubble parameter, and 


oy denotes the error of observed data. 


C. Baryon Acoustic Oscillation (BAO) 


Baryon Acoustic Oscillation provides an efficient method for measuring the expansion 
history of the universe by using features in the cluster of galaxies with large scale survey. 
Here we use the results from the following five BAO surveys: the 6dF Galaxy Survey, SDSS 
DR7, SDSS DR9, WiggleZ measurements and the radial BAO measurement. 
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1. 6dF Galaxy Survey 


The 6dFGS BAO detection can constrain the distance-redshift relation at z,;; = 


0.106 . it gives a measurement of the distance ratio 


faga 
———————  — 0.336 + 0.015 17 
Dy(z = 0.106) ' (17) 
where r,(za) is the comoving sound horizon at the baryon drag epoch when baryons became 


dynamically decoupled from photons. The redshift zg is well approximated by |35 


1291 (Qmoh?)0251 
a rogaz ll 5h)? m 
2 T+ 0.659(Q,.oh®) 055 + bi (Cg h^) ^], (18) 


where 


b, = 0.313(Qmoh 9 49D + 0.607(0,,557)9 674]. (19) 
by = 0.238(Q,,9h2)9 933. 
The effective “volume” distance Dy is a combination of the angular-diameter distance D 4(z) 


and the Hubble parameter H(z), 


1 2 1/3 
Dv) = [t | on 
= [a € 2? DAY gr] 
The xyr is defined by 
a. [(rs(za)/Dy(0.106))+, — 0.336]? 
X6dF — (——— 0018 ^ (21) 


2. SDSS DR7 


The joint analysis of the 2-degree Field Galaxy Redshift Survey data and the Sloan Digital 
Sky Survey Data Release 7 data gives the distance ratio at z = 0.2 and z = 0.35 [30]: 


TRU) ^ 0.1905 + 0.0061, 
Dy(z — 0.2) (22) 
rs(zq) 
— sid) ^ 0.1097 + 0.0036. 
Dy(z — 0.35) 


By applying the reconstruction technique to the clustering of galaxies from the SDSS 
DR7 Luminous Red Galaxies sample, and sharpening the BAO feature, Padmanabhan et 
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al. obtained the distance ratio at z — 0.35 
= 0.1126 + 0.0022. (23) 


The SDSS DR7 and SDSS DR7 reanalysis results are based on the same survey and the 
latter gives a higher precision than the former, we therefore take the SDSS DR7 reanalysis 
data instead of the first one. The x55; ,, used in the Markov Chain Monte Carlo analysis 


is 


2 [Em — 0.1126]? 
Xbre = "9.00007 7 (24) 


3. SDSS DRI 


The SDSS DR9 measurement gives the distance ratio at z = 0.57 "m 


——— ——— —— = 0.0732 + 0.0012. 25 
Dy(z = 0.57) (25) 
The chisquare here is defined as 
rs(Zq) 
2 (DoE) tn — 0.0732]? (26) 
ADRES 0.00122 


4. The WiggleZ measurements 


The WiggleZ team measures the acoustic parameter by encoding some shape information 
on the power spectrum 


(a (27) 


The baryon acoustic peaks measured at redshifts z = 0.44, 0.6 and 0.73 in the galaxy 
correlation function of the final dataset of the WiggleZ Dark Energy Survey give 


A(z = 0.44) = 0.474 + 0.034, 
A(z = 0.60) = 0.442 + 0.020, (28) 
A(z = 0.73) = 0.424 + 0.021. 


The corresponding chisquare is defined as 


Xivig = XTy-x, (29) 
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where 
A(z = 0.44), — 0.474 


X = | A(z = 0.60);, — 0.442] , (30) 
A(z = 0.73), — 0.424 
and its inverse covariance matrix is 
1040.3 —807.5 336.8 
V= = |—807.5 37203 —15519]|. (31) 
336.8 —1551.9 2914.9 


5. Radial BAO 


The radial (line-of-sight) baryon acoustic scale can also be measured by using the SDSS 
data. It is independent from the BAO measurements described above.The measured quantity 
is 

Py (2) em HC ar). (32) 
whose values are given by as 


A, (0.24) = 0.0407 + 0.0011, 


(33) 
A, (0.43) = 0.0442 + 0.0015. 
D. Reduced CMB data 

The chisquare for the reduced CMB data is defined by 

Xomp = X C X, (34) 
where 
(La)th = (La) obs 

X= Rin = Rovs (35) 


(Zath — [Ese onion 


and C is the related covariance matrix. 
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